Cooperative atomic scattering of light from a laser with a colored noise spectrum 
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The collective atomic recoil lasing is studied for an ultra-cold and collisionless atomic gas in a partially 
coherent pump with a colored noise. Compared to white noise, correlations in colored noise are found to be 
able to greatly enhance or suppress the growth rate, above or below a critical detuning. Effects on cooperative 
scattering of light for noise correlation time, noise intensity and pump-probe detuning are discussed. This result 
is consistent with our simulation and linear analysis about the evolution equations in the regions of instability. 
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I. INTRODUCTION 

Collective nonlinear interactions between cold atoms and 
light have attracted considerable attention since collective 
atomic recoil lasing (CARL) was observed in cold atoms |[lll2l 
and superradiance from Bose-Einstein Condensation (BEC) 
was realized in experiments IJ]. The signature of CARL is 
sudden buildup in the probe field or back-scattering beam ori- 
ented reversely to a coherent pump strongly interacting with 
an atomic gas fllUS]. This self-organization phenomenon 
shows spontaneous formation of an atomic density grating, 
and offers the possibility to study the light amplification de- 
rived from collective interaction of light with cold atomic 
gases iHiilii. 

Most studies so far have focused on the coherent laser 
pump, neglecting the fact that environment is intrinsically 
noisy. As found in a recent work for a partially coherent 
pump field with a white noise spectrum jllfl . the probe in- 
tensity can be larger than that in a coherent pump due to 
the noise. The intensity fluctuation and noise play impor- 
tant roles in this nonlocally-coupled many-body system. Al- 
though the white-noise assumption is convenient for mathe- 
matical treatment, it is somewhat unrealistic, because fluc- 
tuations in the microscopic dynamics have a finite (nonzero) 
coiTelation time, commonly referred to as colored noise lfl2ll . 
The white noise approximation is valid when the correlation 
time of fluctuation is much shorter than all other relevant time 
scales in this problem. In other cases there are discrepancies 
between white-noise theory and experiments, such as the pho- 
ton statistics of a dye laser output where the relative intensity 
fluctuations tend to increase indefinitely as the laser is weakly 
excited lfT2l [l3i [3, the reversed asymmetry in the doublet 
spectrum of double optical resonance is reverted back to nor- 
mal for lager detuning ifTsll . Hence, the colored noise model is 
more realistic than the white-noise model and is widely stud- 
ied for laser systems, for example in the sideband squeezing 
in intracavity second-harmonic generation lfT6ll . excess quan- 
tum noise in a laser 1 17i|, and four-wave mixing ifisll . In these 
systems, the colored noise is caused by the laser frequency 
fluctuation, and the field line shape is Gaussian in some pa- 
rameter regime, different with the Lorentzian line shape for 



white noise. 

In this paper we study CARL based on a laser pump with 
colored noise. Eff'ects of noise on the cooperative scattering 
of light are studied by stochastic simulation and linear anal- 
ysis at small correlation time. The Lorentzian bandwidth is 
related to noise intensity, pump-probe detuning, and the cut- 
off frequency given by the inverse of correlation time. Effects 
of these parameters on the growth rate of amplification are 
analyzed and discussed. 



n. MODEL 

We consider the normal setup for CARL The 
backscattered light named probe beam has amplitude Ei and 
frequency oji, and a partially coherent strong pump field has 
amplitude E2 and the mean frequency cl>2- These two beams 
with almost the same frequency |wi - W2I or cl>2, are ap- 

proximately counterpropagating to form a spatially periodic 
optical lattice potential, interacting with a very cold and colli- 
sionless atomic gas. Atoms moving in this lattice form a den- 
sity grating, and in turn play a role on the evolution of probe 
field. 

Since |w2 - <j>iI is assumed to be very small, the pump or 
probe wave vector is given by A: = In /A where A s A2 ~ 
Ai. The pump frequency w s is detuned from the atomic 
resonance frequency wq by A - o) - ajQ. The two-photon 
recoil frequency is given by ojr - 2M~/m. In terms of the 
rescaled time variable r = ojrpt, the rescaled amplitude of 

probe field A - -yJ-^^Ei , the rescaled position 0j - 2kz j and 

momentum pj - {mv j) I {fikp) of the jth atom, the modified 
classical CARL equations with a stochastic pump phase noise 
are given by [1 1] 



dOjldT 
dpj/dr 

dA/dT 

d4>/dT 
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where the dimensionless CARL parameter is given by p = 
(^1^)2/3^ the pump Rabi frequency is Q = dE2/fi, the atomic 
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density is «, the atom-mode coupling constant is given by 
g — d and the dipole matrix element for the atomic 



2 



transition is d. The average in Eq.© is defined as the average 
over all atoms, 

((■■•»-^x|;(...).. 

and the scaled pump-probe detuning is given hy 5 - (0)2 - 
a)\)l{a)rp). 

The partial coherence of pump field is described by a 
phase diffusion model, assumed to evolve according to Eq.©, 
where s {t) i s a Gaussian random variable with zero mean 
<e(r)) s iCr) = 0, and variance |l2l[l8l|20|l 

{s{t)s{t')) = s(T)e(T') ^-e-"^. (5) 

To 

The colored laser noise given by Eq.(5) is parameterized by 
the noise correlation time tq and the noise strength F, where 
(■ ■ ■ ) denotes the average over sampling. In the limit tq — > 0, 
from Eq. (5) we get {^{t)^{t')} — * IFdir - t'), corresponding 
to the case of a white noise in laser phase fluctuations with a 
Lorentzian line shape and Unewidth [half width at half max- 
imum (HWHM)] r for the pump field fTl. When rn -» 00, 
it recovers the case with a coherent pump laser iH [Till . For 
1/to » F, the field line shape is a Lorentzian with a full 
width at half maximun (FWHM) of 2F, whi le for 1/tq <s: F 
it is a Gaussian with a FWHM related to y/r/ro. The gen- 
eral values of Ftq between these two limits lead to Voigt 
profiles iH, meaning that the laser spectrum is significantly 
non-Lorentzian for detuning greater than 1 /tq (from line cen- 
ter) and the spectrum cuts off much more rapidly than a 
Lorentzian. Consequently I/tq is known as the cut-off ifTsll . 
Eq.(5) can also be written as e(t) + -^^(7) - V2r^(T), 
where ^(t) is the Gaussian white noise satisfying (?(t)) = 0, 
{g{T)g(T')} = 6(t - t') iH]. As pointed out in the Ref. ifTTIl . 
the amplification due to the partially coherent pump is less 

sensitive to the atomic distribution in the form of exp(-:^) 
than that due to a coherent pump, where the scaled initial mo- 
mentum is pq, cr — Kg is Boltzmann's constant, and 
T is the temperature of atomic gas. In this paper we consider 
effects of correlation time while assuming cr -Q. 

III. STOCHASTIC SIMULATION 

We directly solve the set of stochastic ordinary differen- 
tial equations ([T]i-(|4li using stochastic Runge-Kutta algorithms 
with a colored noise ll20ll . To conveniently compare with the 
case of a white noise, in the following we use the same value 
F = 5 and detuning i5 = 2,5, 10, as in Ref. 1 1 li] . The scaled 
probe intensity over time is plotted in Fig [l] to show effects 
of pump phase diffusion for pump-probe detuning 6-2. 
The solid line 1 is for the case with white noise F = 5 and 
T"o = ifTIIl . The case of coherent pump with tq = and 
F = is shown in the solid line 2, where the intensity of 
probe beam is very low and oscillates with time, and the gain 
is suppressed. The growth rate and intensity of the backscat- 
tered field greatly increase under the partially coherent pump. 




FIG. 1: (Color online) The scaled probe intensity |Ap (averaged over 
100 runs) for a cold gas (cr = 0) by a partially coherent pump with 
detuning 5 = 2 for different correlation time tq. The solid line 2 is 
for noise intensity F = 0, while all other lines are with F = 5. 
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FIG. 2: (Color online) The scaled probe intensity |Ap versus time for 
detuning (5 = 5,F = 5. 



In the case of colored noise, the growth rates and amplitudes 
of the back scattering beams are increased compared to the 
case of white noise, and this gain increases with correlation 
time To, as shown in the dotted line (tq - 0.02), dashed 
line (to - 0.05), dash-dotted Une (to - 0.08), sohd Une 4 
(to = 0.5), and solid line 3 (to - 2) for 5 = 2. The correlation 
time of noise has an opposite effect with 5 = 5, as shown in 
Fig|2] The growth rates become slow, and it take more time to 
get to the saturation than in the case of 5 = 2. The growth rate 
and intensity of probe field is suppressed for To = 2, as shown 
in the solid line 2. For small correlation time to - 0.02(dotted 
line), it coincides with the case of white noise (the solid line 1 
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FIG. 3: (Color online) The scaled probe intensity |Ap versus time for 
detuning i5 = 10,r = 5. 



To - 0). The growth rate decreases with the increase in corre- 
lation time, as shown in the dashed line (tq = 0.05), and dash- 
dotted line (to = 0.08). However, the difference between the 
colored noise and white noise at these small correlation times 
is not so obvious for 6 = 5. For to = 0.5, there is a distinct 
decrease in the growth rate and intensity of probe beam com- 
paring to the case of white noise, as shown in the solid line 3. 
When the detuning is further increased to ^ = 10, as shown in 
Fig|3] the rate and intensity of gain are greatly suppressed for 
correlation time to = 2 (solid line 2) and tq = 0.5 (solid line 
3). At small To, growth rates also decrease with the increase 
in correlation time, as shown by the dotted line (to = 0.02), 
dashed line (tq = 0.05) and dash-dotted line (tq = 0.08). 

The above results show that the gain is suppressed with the 
increase in correlation time for large detuning (5 = 5 and 
6 = 10). The condition of gain is destroyed with the long 
correlation time of noise. However for the case with small de- 
tuning 5 = 2, the long correlation time helps to enhance the 
amplification in the growth rate and amplitude of the probe 
beam. According to our simulation, there are similar results 
for r = 3 and F = 1 , where a long correlation time is able to 
suppress the growth rate for 6 = 5 and 5=10, and enhance it 
at 5 = 2. Thus depending on detuning, the increase in correla- 
tion time can either enhance or suppress the gain. To further 
understand these diff'erent behaviors, we carry out the linear 
analysis at small correlation time in the next section. 



IV. LINEAR ANALYSIS 



To analyze eff^ects of colored noise in the regions of insta- 
bility, we derive the linear evolution equations for the average 



scattered-field intensity A* A from equations Eqs.([T]i-@, 



dT 

dA*b 

dT 

d'¥P 
dT 

dW 
dT 

db*P 
dT 

dT 



A*b + c.c., (6) 

= -i^ + Vb-iie + T')^, (7) 

= -\A^ + 'FP-(i6 + T')¥P, (8) 

= -ib*P + c.c., (9) 

= -Ab* + i\Pf-, (10) 

= -Fp + c.c, (11) 



where |A|2 = A*A, \b\^ = b*b, \P\- = P*P, and 



N N 



To understand how to derive the above formula, here we 
show the steps in getting Eq.(|7]i. Other equations can be ob- 
tained with similar steps. From Eq.([3]l, we have 

dA*b 



dT 



-iA*P + b*b-i6A*b + isA*b, (13) 



where the last term can be expressed as 



isA*b = /— e-'^^e'P ''(^')''^'/(ti , T), (14) 



with 



/(Tl,r)s<e(T)e^''^<"'"'^'>. (15) 

We can expand Eq.(fT5]l with the series 

^ /(.(tXtOMt, 

( r dT,, f dT,/{e{T,MT,M'- 



X ( 



For small correlation time to, it can be further expressed for 
the moments m = and m = 1 only. 



f(TuT) « iT(\ - e -o)[l + ]. 



Assuming T = t\ + Ktq where K is an adjustable parameter, 
1 «: A" <K (r - Ti), we have 



with 



ieA*b = -r'A*b 



r' = r[i-(/:-i)rTo/2]. 



(16) 



(17) 
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FIG. 4: The growth rate of probe beam %e(A) versus detuning S for 
different correlation times. The solid line 1 is for the coherent pump 
with F = 0, To = 0, and the solid line 2 is for the partially coherent 
pump with white noise F = 5, Tq =0. Other lines are for the partially 
coherent pump with colored noise intensity F = 5 and K = 4, and 
with correlation time Tq = 0.08(dashed line), Tq = 0.05 (dotted line), 
or To = 0.02 (dash dotted line ). 



Eqs. (l6lfTTI) are identical with equations (5)-(10) in 
Ref. ifTHl except that F for the white noise is now replaced by 
r' for the colored noise. Thus we can use the relation between 
the white noise and colored noise models given by Eq.( 17) for 
small T(). An increase of tq can be regarded as decreasing the 
noise intensity or the linewidth of white noise. Hence, for the 
initially cold and collisionless atomic gas, we obtain the char- 
acteristic root A from the solution of Eqs. dSlfTTI). The region 
of instability or amplification satisfies %e{A) > 0. Based on 
the solution, the growth rates of the probe beams %e{A) as 
the detuning 6 for the different correlation time are plotted in 
FigS 

As shown in Fig HI the maximum growth rate occurs at 
^ = 0. In the case of a coherent pump shown by the solid line 
1, the cutoff is so sharp that no instability or amplification of 
the probe beam occurs above 6c - (27/4)'^-' ^ 1.9. The noise 
can extend regions of amplification because the pump phase 
diffusion broadens the region of instability, as shown by the 
solid line 2 for the partially coherent pump with a white noise 
r = 5, To = 0. The case with small tq = 0.02 is very close to 
the case of white noise ifTTIl . For detuning 6 = 2, the growth 
rates increase with increasing tq, as shown in Fig. 4. However, 
as the detuning increases, the growth rates decrease with the 
increase in correlation time for 6 = 5. Compared to the case 
of white noise, for example, the growth rate is increased be- 
low the critical value 6o ~ 3.3, and decreased above it with the 
correlation time tq = 0.05, as shown in Fig.4. Hence, depend- 
ing on detuning, the effect on the growing rate is different for 
different correlation times under the same noise intensity. For 
the different noise strength, the growth rates versus detuning 
with T() = 0.05 are drawn in Fig. 5. We got <5o ~ 1.8 for F = 1, 
^0 ~ 2.7 for F = 3. The critical detuning 6o decreases with 
the decrease in noise strength F. The growth rates are almost 
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FIG. 5: The growth rate of the probe beam 2le(/i) versus detuning 
S for different noise strength. The solid line I is for the coherent 
pump with F = 0, To = 0. The solid line 2 is for white noise with 
F = 1, To = 0, and dotted line is for colore noise with F = 1, Tq = 
0.05. The solid line 3 is for F = 3, To = 0, and dashed line is for 
F = 3,To = 0.05. 



same for white noise and colored noise with small noise in- 
tensity F = 1, and there is an obvious difference for F = 3. 
The increase in the correlation time is beneficial to the growth 
rate of the amplification near the threshold detuning, while 
destructive for the big detuning. 



V. DISCUSSION AND CONCLUSIONS 

The gain of CARL given by the growth rate of the probe 
beam depends on the frequency of the light field. It has a cer- 
tain gain bandwidth defined as the spectral range where the 
light scattering is exponentially amplified. In the amplifica- 
tion process, atom scatter photons from the pump field into 
the probe beam, experience an acceleration due to the pho- 
tonic recoil, and occupy different momentum states. On the 
other hand, the frequency of scattered photons is shifted with 
respect to the pump-light frequency, resulting in pump-probe 
detuning. In the case of a noisy pump, there are two velocities 
involved. One is the stochastic phase velocity with mean value 
{a)[ - a)2)c/{a>i + 0)2) of the optical potential, and the other is 
the atomic velocity. Due to the interaction between the atoms 
and the optical field, the synchronization between these two 
velocities eventually leads to the scaled momenta distribution 
of atoms around a mean value {p} ^ -6, when 6 > 6c = 1.9 
and the dynamical phase evolution of probe field is negligi- 
ble. Thus the region of instability extends from the threshold 
value of the coherent pump 6c to the bigger detuning ifTHl . At 
the same time, a decrease in (p) causes an increase in |Ap, due 
to the conservation of momentum (p) -H |Ap = const. Hence 
the saturation of instability occurs when the scaled probe in- 
tensity satisfies |Ap » 6. This is the reason of amplification in 
large detuning for the partially coherent pump. 

However, the noise correlation time tq, noise intensity F, 
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and pump-probe detuning 6 > 6c greatly affect the growth rate 
of this ampHfication behavior. To understand this relation, we 
define y as the ratio of the effective gain bandwidth Te/f deter- 
mined by noise intensity and correlation time to the detuning, 

r = Tcff/S. (18) 

It determines how many momentum states lying within the 
effective linewidth may participate and be amplified in the 
CARL dynamics. If this ratio is one, it means that all the mo- 
mentum states within the linewidth are amplified at the same 
time and the growth rate is the fastest. If the ratio is more or 
less than one, this growth rate decreases. 

In the case with small correlation times, F <s; I/tq, the 
cut-off of the laser spectrum is larger than Lorentzian band- 
width. The colored noise model and white noise model lead 
to similar results following the relation given by Eq.(17), 
Teff - r[l - (K - l)rTo/2]. The effective bandwidth of noise 
Feff limits the range of frequencies accessible for the probe 
light field. For r = 5 and 5 = 2, r a! 2.5 - 18.75to, the 
increase in correlation time leads to the decrease in y closer 
to one. Hence the growth rate of amplification is enhanced 
with increasing tq, as shown in Fig. 4 of linear analysis and in 
Fig.l of simulation results with tq - 0.02,0.05,0.08. How- 
ever, for S - 5, y ~ 1 - 7.5to, the increase in correlation 
time leads to a deviation from the match condition, hence in- 
creasing To decreases the growth rate, as shown in Fig.4 of 
the linear analysis and in Fig. 2 of simulation results. Further- 
more, for (5 = 10, y a; 0.5 - 3.75to, increasing correlation 
time destroys the match condition y - 1, as demonstrated in 
Fig.4 of linear analysis and in Fig. 3 of the simulation results 
with To = 0.02, 0.05, 0.08 . Based on linear analysis shown in 
Fig.4, the simulation results at small correlation times shown 
in Fig. (1-3), and the above analysis, we reach the conclu- 
sion that the change in the correlation time can be regarded 
as adjusting the effective noise intensity. Whether the change 
in correlation time enhances or suppresses the growing rate 
depends on whether it helps or destroys the match condition 

r= 1. 

The large correlation time corresponds to slow frequency 
fluctuation I/tq <*c F, when the laser shape approaches a 
Gaussian with the effective bandwidth of HWHM Fg/y - 
V[8(/n2)r/To]/2. The effective bandwidth limits the range of 
frequencies accessible for probe beam. When 5 = 2, we have 
y - 1.86,0.93 for tq - 0.5 and 2. Because 0.93 is closer to 
unity than 1.86, the growth rate at tq = 2 is faster than that at 
Tq = 0.5, as shown in the solid hne 3 and 4 of Fig. 1 . However, 
for 5 = 5, To = 0.5, 2 corresponding to y = 0.74, 0.37, the in- 
crease of correlation time destroys the match condition and the 
amplification at to = 2 is almost suppressed, as shown in the 
sold line 2 of Fig. 2. Furthermore, at 5 - 10 and to = 0.5,2, 
y = 0.37, 0.18, these two values are far from one and the am- 
plification is also suppressed at tq = 0.5, as shown in the solid 
line 3 of Fig. 3. The simulation results show that the growth 
rate is greatly suppressed with large correlation times. Hence, 
for big To, the increase in correlation time greatly reduces the 
bandwidth, and suppresses the number of momentum states of 
participating in the CARL dynamics. The above explanation 



also is suitable for F = 3, while it is not right for weak noise 
intensity F = 1 . 

It should be possible to observe these effects following 
the experimental observation of CARL iQ], |2] where ultra- 
cold ^^Rb atoms were enclosed in a ring cavity. The differ- 
ence between the experimental setup and our model is that 
the pump field in the experiment is also in a cavity mode 
which counterprop agates with respect to the probe beam. 
The experimental reported characteristic growth time tg for 
the instability of CARL is about l^s corresponding to 
a>,.p a; f-' = 10*'i"', p a 10, and the recoil shift of Rb 
as a>r = 2M~/m ^ In x lAkHz- For a high finesse of the 
cavity 6u)pump ^ In x IQkHz and the scaled pump linewidth 
F = 6cL>p,„„p/(a)rp) ~ 0.02, it corresponds to the case of a 
coherent pump field for F «: 1 . However, the linewidth of 
the pump field can be adjusted between good-cavity and bad- 
cavity regimes by varying the finesse of cavity, atom num- 
ber, and pump power. F from 0.3 to 0.2 corresponding to 
p - 4.7 to 7.0 is the semiclassical good-cavity regime. F 
from 3.7 to 2.8 corresponding to the parameter p - 5.1 to 
6.7 is the typical semiclassical bad-cavity regime IH 01- If 
we assume p = 5 for Rb, the scaled noise intensity F = 5 
means a linewidth of the Lorentzian of In x 350kHz, the 
scaled correlation time to - 0.5 means the cutoff at frequency 
Itt X 14-QkHz, and pump-probe detuning 6 - 2,5, 10 corre- 
sponds to 14-OkHz, 35QkHz, lOQkHz, respectively. Hence, to 
certain extent, our analysis can be regarded as the CARL in 
bad-cavity regime which is possible to be accessed in experi- 
ments. 

In conclusion, the correlation time of pump phase noise 
greatly affects the growth rate and intensity of cooperative 
scattering in the system of CARL. The noise makes the am- 
plification region for pump-probe detuning larger than that in 
the coherent pump. The noise intensity and correlation time 
determine the effective linewidth of pump laser, not just noise 
intensity in the case of white noise predicted by phase dif- 
fusion model of ideal laser theory. The change in correla- 
tion time can enhance or suppress the growing rate depending 
on the ratio of the effective bandwith to pump-probe detun- 
ing which determines how many momentum states within the 
linewidth are amplified at the same time. This ratio equal to 
one corresponds to the best match condition for big noise in- 
tensity. Whether the growth rate is suppressed or enhanced by 
the change in correlation time depends on whether it helps or 
destroy the synchronization condition. These results are use- 
ful for analyzing the cooperative scattering process and effects 
of noise on the collective nonlinear interaction between cold 
matter and fight, and also helpful in studying superradiance 
from BEC lulll or the phase coherent matter-wave amplifica- 
tion ll22ll because of the same gain mechanism. 
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